GoS i 


Time Dependent Density Functional Theory Calculations of Large 
Compact PAH Cations: Implications for the Diffuse Interstellar 

Bands 

Jennifer L. Weisman 1 , Timothy J. Lee 2 , Farid Salama 3 and Martin Head-Gordon 1 ' 4 


ABSTRACT 


We investigate the electronic absorption spectra of several maximally peri- 
condensed polycyclic aromatic hydrocarbon radical cations with time dependent 
density functional theory calculations. We find interesting trends in the verti- 
cal excitation energies and oscillator strengths for this series containing pyrene 
through circumcoronene, the largest species containing more than 50 carbon 
atoms. We discuss the implications of these new results for the size and structure 
distribution of the diffuse interstellar band carriers. 

Subject headings: astrochemistry — ISM: molecules methods: analytical 


1. INTRODUCTION 

Polycyclic aromatic hydrocarbons (PAHs) are commonly accepted as being ubiquitous 
in the interstellar medium (ISM). They have received considerable attention in recent years 
as possible carriers of the diffuse interstellar bands (DIBs), after the presentation of the 
PAH hypothesis (Crawford et al. 1985; van der Zwet and Allamandola 1985, Leger and 
d’Hendecourt 1985; Salama et al. 1996, 1999). There have been many studies both experi- 
mental and theoretical, on various types of PAHs that may be DIB carriers, and that may be 
most likely to exist in the ISM (for some example references, see Salama et al. (1994, 1995); 
Joblin et al. (1999); Romanini et al. (1999); Brechignac and Pino (1999); Weisman et al. 
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(2001)). PAH cations, both radical and closed-shell, are good candidates for DIB carriers. 
Not only do they have electronic absorptions in the spectral region of interest, but they ad- 
ditionally have high ionization energies, a requirement to survive in the ISM (Sonnentrucker 
et al. 1997). 

Theoretical predictions of PAH absorption spectra are a useful tool in the search for DIB 
carriers. Ideally, gas-phase spectroscopy in the laboratory is necessary to be able to match 
accurately a potential DIB carrier to astrophysical data. However, current experimental dif- 
ficulties in sample preparation have inhibited direct gas phase spectroscopy of larger PAHs, 
and most studies have been performed in rare gas matrices. Although recent studies (Ro- 
manini et al. 1999; Brechigr.ac and Pino 1999; Biennier et al. submitted 2002) have provided 
the first electronic spectra of PAH ions in the gas phase, to our knowledge, absorption results 
have not yet been published for medium to large sized PAHs in the gas-phase (i.e. PAHs 
containing more than 20 carbon atoms). In this regard, calculated spectra are potentially 
useful in guiding future experimental work. Calculated spectra are also valuable in terms of 
gai ning a better understanding of the nature of the transitions (Gutman and Cyvin 1989). 

Time dependent density functional theory (Runge and Gross 1984; Petersilka et al. 
1996) methods are capable of predicting vertical excitation energies and oscillator strengths 
for relatively large molecules (for example, see Weisman et al. (2001)). Excitation energy 
accuracy is possible with conventional exchange-correlation functionals to within less than 
a third of an electron volt (eV) for radical cations such as pyrene (Hirata et al. 1999). The 
ability to predict oscillator strengths is also quite useful. Our experience is that oscillator 
strength predictions are in qualitative accord with experimental observation (Hirata et al. 
1999). 

In a recent study, we have found an interesting trend in oscillator strengths within the 
perylene, terrylene, quaterrylene series (Halasinski et al. submitted 2002). The bright state 
for each of these molecules is a transition of dipole symmetry along the long axis of the 
molecule. Since the dipole forms along the long axis, the oscillator strength of the bright 
state increases as the molecular backbone grows in length from perylene to terrylene to 
quaterrylene. This results in a net increase in oscillator strength per unit mass of carbon 
along this series. Such a conclusion could not be easily made without TDDFT calculations, 
since an oscillator strength is very difficult to quantify in the only experimental techniques 
that have so far been used o study these species, namely matrix-isolation spectroscopy. 

This information has potentially significant astrophysical implications. A larger oscil- 
lator strength per unit mass of carbon for larger molecules could certainly shift the search 
for DIB carriers towards larger species if this trend holds for PAHs other than just this 
one oligorylene series. This is what we aim to find out by studying the oscillator strength 
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trends in more PAH series In this letter, we have chosen to investigate a series of maxi- 
mally pericondensed PAH radical cations. Maximally pericondensed PAHs are those PAHs 
that have the largest number of carbon atoms for a given number of hydrogen atoms. These 
PAHs have received much attention since they are expected to be especially stable structures 
(Stein 1978) and the likely precursors to soot formation in hydrocarbon flames (Frenklach 
and Feigelson 1989; Cherchneff and Barker 1992). 


2. ELECTRONIC STRUCTURE METHODS 

All calculations were performed with the Q-Chem quantum chemistry program package 
(Kong et al. 2000) on IBM RS/6000 workstations. The following PAH radical cations were 
investigated: pyrene (Ci 6 H 0 ), coronene (C 24 H 12 ), ovalene (C 32 H 14 ), circumpyrene (C 42 Hi 6 ), 
and circumcoronene (C 54 H ; g). See Figure 1 for oriented structures. All ground state geome- 
tries were optimized using density functional theory with Becke exchange (Becke 1988) and 
LYP correlation (Lee et al. 1988), and the 6-31G* basis set. The geometries were confirmed 
to be minima by calculation of the Hessian matrix. 

The vertical excitation energies and oscillator strengths were determined using time 
dependent density functional theory (TDDFT) at the BLYP/6-31G* level. Radical cation 
excitation energies are generally insensitive to basis set size for TDDFT calculations (for 
example, see Halasinski et al. (submitted 2002)), so we believe the 6-31G* basis is sufficient. 
Commonly used LDA and GGA functionals perform well for predicting the energies of va- 
lence states, but perform quite poorly for Rydberg states (Casida et al. 1998; Tozer and 
Handy 1998; Hirata and Head-Gordon 1999). However, we roughly estimate the onset of 
the Rydberg level as the negative of the highest occupied Kohn-Sham orbital (Casida et al. 
1998), which falls in the range of 7-9 eV for the PAH radical cations in this study. Our 
reported excitation energies remain well below this level. 


3 RESULTS AND DISCUSSION 

The PAH radical cation calculations in this investigation have all been run with D 2h 
molecular symmetry. Neutral coronene and circumcoronene have D 6/l molecular symmetry, 
however the radical cation derivatives are expected to experience Jahn-Teller effects and 
distort to D 2h symmetry (Kato et al. 1999). We orient the pyrene radical cation in the xy 
plane with the long axis along the x axis. All other molecules contain a central ’pyrene- 
like’ unit with the same or entation, so that for larger molecules like ovalene, the long axis 
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is actually in the y direction. The coronene and circumcoronene D 2 h structures are very 
slightly elongated in the y direction. 

We focus on the brightest transitions of the series, that is, the electronic excitations with 
the largest oscillator strengths. These are the transitions with most promise for correlation 
with the DIBs (Salama et ai. 1996, 1999). We present the strongest excitation of B 2u and B 3u 
transition symmetry. The transition symmetry is the product of the ground state and excited 
state symmetries, and is a simpler way to characterize the transitions m these molecules, 
since the ground state symmetries are different. In the xy plane orientation described earlier, 
the B 2u and B 3u transition? are oriented along the y and x axes, respectively. 

Table 1 shows the bright transitions for the pyrene, coronene, ovalene, circumpyrene and 
circumcoronene radical cations. It is interesting to note that coronene and circumcoronene 
do not seem to fit into a trend with the other radical cations. The B 2u excitation is lower 
in energy than the B 3u excitation for all of the molecules. However, the transition with the 
larger oscillator strength is not the same for all five species. B 3tl is the stronger absorption 
for pyrene, ovalene and circumpyrene, while B 2u is the stronger absorption for coronene and 
circumcoronene . 

These new theoretical iesults confirm the previous findings regarding the general spectral 
characteristics of PAH ions (Salama et al. (1996), and references therein). They show that the 
visible absorption spectra of pericondensed PAH ions are indeed dominated by one transition 
although the relative intensities can differ among different subfamilies. They also provide, 
for the first time, precise information about the nature of the molecular transitions as well as 
important input into the quantitative variation of the oscillator strengths with the molecular 
structure of the PAH ions. Oscillator strengths are key parameters for astrophysical models 
because they dictate the abundance of cosmic carbon that needs to be locked-up in any 
specific molecular carrier, and, hence, its plausibility as an interstellar band carrier. Armed 
with these new results, one can now fine-tune the theoretical predictions that were made in 
previous astrophysical models (Salama et al. 1996). 

It is clear that circumpyrene behaves similarly to pyrene, and likewise circumcoronene 
mirrors coronene. Let us refer to the pyrene-circumpyrene pair as the pyrene series, and 
the coronene-circumcoroneae pair as the coronene series. The next larger member of these 
homologous series is then defined by adding benzene rings to completely enclose circumpyrene 
and circumcoronene, analogous to the way that circumpyrene and circumcoronene were 
formed from pyrene and coronene, respectively. 

In both the pyrene and coronene series, the larger species has an excitation energy 
roughly 0.6 eV lower than ,he smaller species for the B 2u transition, and 1 eV lower for the 
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B 3u transition. Additionally, the oscillator strength of the larger species in relation to the 
smaller species increases for the B 2 „ transition and decreases for the B 3u transition. The 
pyrene series strong absorption (B 3u ) becomes weaker for larger members, so the oscillator 
strength per unit mass of carbon quickly decreases. The B 2u absorption does grow stronger 
for larger members, however it is fairly weak to begin with and the excitation energy itself 
decreases rapidly out of the range of the known DIBs. The situation is similar for the 
coronene homologous series, except that the B 3u absorption is not strong. We will discuss 
these trends further later. 

Ovalene behaves similarly to pyrene and circumpyrene. In comparison to circumpyrene, 
ovalene’s oscillator strength increases more relative to pyrene for the B 2u transition and 
decreases more from pyrene for the B 3u transition. This is most likely due to the fact that 
the long axes for ovalene and circumpyrene are switched, that is ovalene’s long axis is the 
y axis, and so the B 2u transition is stronger than circumpyrene’s, and visa versa for the 
B 3u transition. Overall, ovalene seems closer to the pyrene series than the coronene series, 
however it is not part of a homologous series in either case so we will not pursue it further 
in this discussion except to emphasize that all compact PAHs may not be grouped together 
for purposes of predicting oscillator strength trends. 

Let us consider the pyrene series and coronene series further now. We have already seen 
that the oscillator strength of the B 2u transition decreases along both homologous series and 
the strength of the B 3u absorption increases. What is even more important to consider in the 
context of the interstellar carbon budget, is the oscillator strength per unit mass of carbon. 
Table 2 shows the oscillator strength and oscillator strength per unit mass of carbon for 
both transitions. We can now see that while the oscillator strength itself increases for the 
B 3u transition, the oscillator strength per unit mass of carbon does not. In other words, the 
species in these series grow in size faster than the strength of their absorptions. Thus the 
oscillator strength per unit mass of carbon decreases for both of the strongest absorptions 
in the pyrene series and the- coronene series. 

The results for these two series of PAH radical cations show trends that would seem 
to argue against larger molecules of these series as promising candidates for DIB carriers, 
despite their stable structures. With the exceptions of coronene and circumcoronene, these 
PAH radical cations do exhibit transitions of strong oscillator strength. However, these 
strong transitions decrease per unit mass of carbon along the series and except for pyrene, 
are not comparable to the perylene series. For example, the oscillator strength per unit 
mass of carbon for quaterrylene’s strongest absorption is approximately eight times larger 
than circumcoronene ’s strongest absorption. This would require almost 10 times more cir- 
cumcoronene in the intersl ellar medium than quaterrylene to yield a signal of the same 
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magnitude. Additionally, the excitation energies of the pyrene and coronene series B 2u tran 
sitions decrease rapidly out of the range of the known DIB region (roughly 1.0 - 3.1 eV). 


4. CONCLUSIONS 

We have presented calculated absorption spectra of several pericondensed PAH radical 
cations in this letter. Time dependent density functional theory methods allow for the 
investigation of larger molecules than are readily feasible with other ab initio electronic 
structure methods. We have studied the first five maximally pericondensed PAH radical 
cations: pyrene, coronene, ovalene, circumpyrene, and circumcoronene. 

We have found interest ing and insightful trends in the oscillator strength per unit mass 
of carbon for two series of pericondensed PAH radical cations. First, we have identified the 
fact that oscillator strength trends must be addressed in specific homologous series only. It is 
not appropriate to generalize trends for all compact or non-compact PAHs together because 
they may behave very differently. We have identified two homologous series within the 
maximally pericondensed PAHs family, namely the pyrene series and coronene series, which 
are grouped with their cir<um- counterparts. These two series differ in that the strongest 
absorption is along the x axis for pyrene (B 3u ) and along the y axis for coronene (B 2u ). 
However, both series exhibit the behavior that the two strongest absorptions in the DIBs 
spectral region decrease in oscillator strength per unit mass of carbon along the series. 

We therefore find the pyrene and coronene radical cation series to be poor candidates for 
carriers of the diffuse interstellar bands in comparison to strong absorbers like the oligorylene 
series, although pyrene itself is still a possible candidate. The decreasing oscillator strength 
per unit mass of carbon and rapidly decreasing excitation energy both contribute unfavorably 
to the possibility of these compact PAHs being DIBs carriers. We should note that we have 
confined our study to the brightest states of roughly the visible region of the spectrum, which 
is the spectral region of interest for the known DIBs. There certainly exist higher lying 
excited states for these PAHs that absorb in the UV and are simply out of the energetic 
range of the known DIBs. 

Our results have specific implications for future work. We have identified that different 
series of compact PAHs behave differently from each other. It seems appropriate that future 
effort should be directed to study the oscillator strengths of a variety of PAH homologous 
series in the hope of identifying both good and bad candidates for the DIBs carriers to aid 
in the direction of future laboratory investigations. We aim to do this in follow up papers. 
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Table 1: Pericondensed Series Brightest Electronic_^VaDsitions oPBg^&nd Bg u SymnietTy 


Radical Cation 

Ba* 

g - 

Transition 

B 3u 

Transition 



Energy 

OS 


Energy 

OS 

Pyrene 

2 B3 g — 

1.5923 

(0.0156) 

2 B3 5 “> 2 A u 

2.9223 

(0.2019) 

Coronene 

2 B lu -» 2 B 3ff 

1.6481 

(0.0473) 

2 B lu -> 2 B 2fl 

2.8391 

(0.0267) 

Ovalene 

2 B 2ff — > A„ 

1.1801 

(0.0293) 

2 B 2fl — > 2 Bi u 

2.1110 

(0.1319) 

Circumpyrene 

2 Bi u — kB 3s 

0.9898 

(0.0263) 

2 B\ U ”> 2 B2 g 

1.9779 

(0.1466) 

Circumcoronene 

2 B 3s — eBi„ 

1.0108 

(0.0660) 

2 B Sg — > 2 A U 

1.8925 

(0.0262) 


a 


BLYP/6-3lG*//BLYP/6-3lG* excitation energies in eV, oscillator strengths (OS) given in parentheses. 
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Table 2: Oscillator Strength (OS) Trends of Pericondensed Transitions of B 2u and B 3u Sym- 


metrv ===^=^= 

Radical Cation 

No. of 

B 2u Transition 

B 3u Transition 


Carbon Atoms 

OS 

OS/Carbon*1000 

OS 

OS/ Carbon* 1000 

Pyrene 

16 

0.0156 

0.98 

0.2019 

12.6 

Circumpyrene 

42 

0.0263 

0.63 

0.1466 

3.49 

Coronene 

21 

0.0473 

1.97 

0.0267 

1.11 

Circumcoronene 

51 

0.0660 

1.22 

0.0262 

0.49 
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Fig. 1. — Oriented Structures of the Pericondensed PAH Radical Cations: 
^Coronene, c) Ovalene, ^Circumpyrene, and e) Circumcoronene. 
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